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Abstract. In this work we study the collective behaviors arising in a 
model of a simplified homogeneous society. Each agent is modeled as a 
binary "perceptron" , receiving neighbors' opinions as inputs and out- 
putting one of two possible opinions, according to the conformist atti- 
tude and to the external pressure of mass media. For a neighborhood 
size greater than three, the system shows a very complex phase diagram, 
including a disordered phase and chaotic behavior. We present analytic 
calculations, mean fields approximation and numeric simulations for dif- 
ferent values of the parameters. 



1 Introduction 

In recent years there has been a great interest in the applications of physical 
approaches to a quantitative description of social and economic processes [1,2]. 
The development of the interdisciplinary field of the "science of complexity" has 
lead to the insight that complex dynamic processes may also result from simple 
interactions, and even social structure formation could be well described within 
a mathematical approach. 

In this paper we are interested in modeling the human society in order to 
understand the process of opinion formation applied to political elections. Many 
different approaches have been taken by scientists in this matter. The simplest 
choice is to consider a uniform society, where all the persons are considered equal 
and have the same influence on all the others. The second choice is to consider a 
society divided into two groups, a group of normal people and a group of leaders 
(see Ref. [3]). 

Another possibility is to consider all the individuals different from one an- 
other. This can be the case of the Kauffman model, where every individual has 
different (randomly generated) features and interacts randomly with a certain 
number of other individuals. Despite its simplicity, this model is able to exhibit 
very complex behaviors, including chaos. 



In this paper we start exploring the behavior of a model of an uniform society. 
By uniform we mean that everybody is exactly like the others, except for the 
expressed opinion (that may vary with time). All other parameters are exactly 
the same. 

We use a simple neural network model to approximate the process of opinion 
formation taking place in a society. With this approximation, the individual 
is represented by an automaton that receives some external inputs, elaborates 
them using a certain function, and elaborates a response. 

As a first step, we study the behavior of the simplest, one-dimensional neural 
network, composed by binary perceptrons. Such a model has been successfully 
applied to anticipate personal preferences on products (see Ref. [4]). The inputs 
for each perceptron are given by the opinions expressed by other persons in a 
local community, whose size is R. 

The parameters of the model are the influence of the mass- media (H), the 
weight J assigned to the the local community (that can be thought as the result of 
education), and the conformist threshold Q. This last parameter models the em- 
pirical fact that even a strong mass-media campaign or a strong anti-conformist 
attitude cannot modify an opinion if it is supported by a strong majority in the 
local community. If the local majority is outside the thresholded intervals, the 
evolution rule is that of a Monte-Carlo simulation of an equilibrium system with 
heat-bath dynamics. 

The system may be trapped into one of the two absorbing states (uniform 
opinion), or exhibit several kinds of irregular behavior. Technically, it is an ex- 
tension of a probabilistic cellular automata with a very rich phase diagram [5], 
whose mean- field approximation was described in Ref. [6] . A detailed investiga- 
tion of its behavior will be illustrated elsewhere [7]. 

In Sec. 2 we present this probabilistic cellular automata model and present 
the main results. In Sec. 3 we extend the concept of Lyapunov exponents to 
probabilistic cellular automata and present the results obtained for our model. 
Finally, in Sec. 4, we summarize our work and draw some conclusions. 

2 The model 

Let xl be the opinion assumed by individual i at time t. As usual in spin mod- 
els, only two different opinions, -1 and 1, are possible. The system is made up 
of L individuals arranged on a onc-dimcnsional lattice with periodic boimdary 
conditions. Time is considered discrete (corresponding, for instance, to election 
events). The state of the whole system at time t is given by a;* = (xq, . . . , a;^_i) 
with xl G {-1,1}. 

The individual opinion is formed according to a local community "pressure" 
and a global influence. In order to avoid a tie, the local community is made up 
of i? = 2r + 1 individuals, including the opinion of the individual himself at 
previous time. The average opinion of the local community around site (person) 
i at time t is denoted by m = m* = . 



Let J be a parameter controlling the influence of the local field in the opinion 
formation process and H be the external social pressure. One could think of H as 
the television influence, and J as the educational effects. The "field" H pushes 
toward one opinion or the other, and people educated toward conformist will 
have J > 0, while non-conformists will have J < 0. 

The hypothesis of alignment with an overwhelming local majority is rep- 
resented by the parameter Q, indicating the critical size of local majority. If 
m<2Q-R, then a;*+^ = -1, and if m > i? - 2Q, then = 1. 

For simpler numerical implementation and plotting, the opinions (—1,1) are 
replaced by (0,1). Let us denote by S the sum of opinions in the local community 
using the new coding. The "magnetization" can by expressed as m = 2S — R, 
and the probability Ps = P(1|S') of expressing opinion 1 given S neighbors with 
opinion 1 is: 



For Q — the model reduces to an Ising spin system. For all Q > we 
have two absorbing homogeneous states, a; = and x = 1, corresponding to 
an infinite plaquctte coupling in the statistical mechanical sense. With these 
assumptions, the model reduces to a one-dimensional, one-component, totalistic 
cellular automaton with two absorbing states. 

The order parameter is the fraction c of people sharing opinion 1^. It is 
zero or one in the two absorbing states, and assumes other values in the active 
phase. The model is symmetric since the two absorbing states have the same 
importance. 

The quantities H and J range from —oo to oo. For easy plotting, we use 
the parameters J = [1 + exp(— 2J)]~^ and h — [1 + exp{—2H)]^^ as control 
parameters, mapping the real axis (— oo, oo) to the interval [0, 1]. 

The fraction c of ones in a configuration and the concentration of clusters p 
are defined by 



In the mean-field approximation, the order parameters c and p are related by 
p = 2c(l — c). Both the uniform zero-state and one-state correspond to p ^ in 
the thermodynamic limit. 

Typical patterns for ff = are shown in Fig. 1. Roughly speaking, for ferro- 
magnetic (J > 0) coupling, the only stable asymptotic states are the absorbing 
ones. The system quickly coalesce into large patches of zeroes and ones (Fig. 1-D), 
whose borders perform a sort of random motion until they annihilate pairwise. 
For J < the stable phase is represented by an active state, with a mixture of 

^ The usual order parameter for magnetic system is the magnetization M = 2c — 1. 
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Fig. 1. Typical patterns (256 x 256 pixels) for RllQl and H = 0, starting from a 
disordered initial condition with po = 0.5. (A) [top left, j = 0.056250]: "chaotic" 
phase. One can observe rare "triangles" with "base" larger that R, corresponding to the 
unstable absorbing states, and other metastable patterns, corresponding to additional 
absorbing states for J — > — oo. (B) [top right, j = 0.421875]: active phase, the only 
absorbing states are and 1. (C) [bottom left, j = 0.478125]: disordered phase. (D) 
[bottom right, j = 0.562500]: quiescent phase. In this phase the only stable states 
arc the absorbing ones. The boundaries separating the phases move randomly until 
coalescence. 
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Fig. 2. Behavior of p for // = 0. RiQl (long dasiied), iillQl (solid), RAIQI (dashed) 
and R&lQl (dotted). Letters correspond to patterns in Fig 1. 



zeroes and ones. For J ^ oo the automaton become deterministic, of "chaotic" 
type (for R=3 it is Wolfram's rule 150). 

As illustrated in Fig. 2, when R is greater than 3, the quantity p is no more 
a monotonic function of j, and a new, less disordered phase appears inside the 
active one for small values of j. This phase is characterized by a large number of 
metastable states, that become truly absorbing only in the limit J ^ — oo. For 
reasons explained in the following, we shall denote the new phase by the term 
irregular, and the remaining portion of the active phase disordered. By further 
increasing R, both transitions become sharper, and the size of the disordered 
phase shrinks, as shown in Fig. 2. By enlarging the transition region one sees 
that for iJ = (Fig. 3) the transitions are composed by two sharp bends, 
which are not finite-size or time effects. As shown in Fig. 1-C, in this range 
of parameters the probability of observing a local absorbing configurations (i.e. 
patches of zeroes or ones) is vanishing. All others local configurations have finite 
probability of originating zeroes or ones in the next time step. The observed 
transitions are essentially equivalent to those of an equilibrium system, that 
in one dimension and for short-range interactions cannot exhibit a true phase 
transition. The bends thus become real salient points only in the limit R ^ oo. 

The origin of the order-disorder (II) and disorder-irregular (I) phase transi- 
tion is due to the loss of stability of the fixed point c* ^ 0,1 given by the mean 
field approximation. More detailed considerations can be found in [7] . 

For H — (Fig. 3) the period-doubling phase brings the local configuration 
into an absorbing state, and the lattice dynamics is therefore driven by inter- 
actions among patches which are locally absorbing. This essentially corresponds 
to the dynamics of a Deterministic Cellular Automata (DCA) of chaotic type. 
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Fig. 3. Comparisons between numerical (tliin line) and mean-field (thick dotted line) 
results for RS>1Q1 and H = {h — 0.5). The estimated critical values are jf ~ 0.493827 
and ju ^ 0.51384. 

i.e. a system which is insensitive of infinitesimal perturbations but reacts in an 
unpredictable way to finite perturbations. 

3 Chaos and Lyapunov exponent in cellular automata 

State variables in cellular automata are discrete, and thus the usual chaoticity 
analysis classifies them as stable systems. The occurrence of disordered patterns 
and their propagation in stable dynamical systems can be classified into two 
main groups: transient chaos and stable chaos. 

Transient chaos is an irregular behavior of finite lifetime characterized by the 
coexistence in the phase space of stable attractors and chaotic non attracting sets 
- namely chaotic saddles or repellers [8]. After a transient irregular behavior, 
the system generally collapses abruptly onto a non-chaotic attractor. Cellular 
automata do not display this kind of behavior, which however may be present 
in systems having a discrete dynamics as a limiting case [9]. 

Stable chaos constitutes a different kind of transient irregular behavior [10, 
11] which cannot be ascribed to the presence of chaotic saddles and therefore to 
divergence of nearby trajectories. Moreover, the lifetime of this transient regime 
may scale exponentially with the system size (supertransients [10, 11]), and the 
final stable attractor is practically never reached for large enough systems. One is 
thus allowed to assume that such transients may be of substantial experimental 
interest and become the only physically relevant states in the thermodynamic 
limit. 

The emergence of this "chaoticity" in DCA dynamics, is effectively illustrated 
by the damage spreading analysis [12,13], which measures the sensitivity to 



initial conditions and for this reason is considered as the natural extension of the 
Lyapunov technique to discrete systems. In this method, indeed, one monitors 
the behavior of the distance between two replicas of the system evolving from 
slightly different initial conditions. The dynamics is considered unstable and the 
DCA is said chaotic, whenever a small initial difference between replicas spreads 
through the whole system. On the contrary, if the initial difference eventually 
freezes or disappears, the DCA is considered non chaotic. 

Due to the limited number of states of the automata however, damage spread- 
ing does not account for the maximal "production of uncertainty" , since the two 
replicas may synchronize locally just by chance (self-annihilation of the dam- 
ago). Moreover, there are different definitions of damage spreading for the same 
rule [14]. 

To better understand the nature of the active phase, and up to what extent 
it can be denoted chaotic, we extend the finite-distance Lyapunov exponent 
definition [15] to probabilistic cellular automata. A similar approach has been 
used in Ref. [16], calculating the Lyapunov exponents of a Kauffman random 
boolean network in the annealed approximation. As shown in this latter paper, 
this computation gives the value of the (classic) Lyapunov exponent obtained 
by the analysis of time-series data using the Wolf algorithm. 

Given a Boolean function f{x, y, . . . ), we define the Boolean derivative df /dx, 

as 

dl^il if/(|x-l|,y,...)^/(x,y,...), 
dx \Q otherwise, 

which represents a measure of sensitivity of a function with respect to its argu- 
ments. The evolution rule of a probabilistic cellular automaton may be thought 
as a Boolean function that depends also by one or more random arguments. In 

our case 

f{xr,X2...;r) = f{S;r) = [r<Psl 

where S = X\ + X2, ■ ■ ■ , Ps '^s, defined by Eq. (1) and [statement] is the truth 
function, which gives 1 if the statement is true and otherwise. In this case the 
derivative is taken with respect to Xi by keeping r constant. 

For a cellular automaton rule, we can thus build the Jacobian Jij = dx*'^^ / i9a;* . 
This Jacobian depends generally on the point in phase-space (the configurations) 
belonging to a given trajectory. In the case of a probabilistic automaton, the tra- 
jectory also depends on the choice of the random numbers r. 

Finally, the maximum Lyapunov exponent A is compiitod in the usual way by 
measuring the expansion rate of a "tangent" vector v{t), whose time evolution 
is given by 

v{t + l) = Jv{t). 

As explained in Ref. [15], a component Vi{t) of this tangent vector may be 
thought as the maximum number of different paths following which any ini- 
tial {t = 0) damage may reach the site i at time t, without self-annihilation. 
When all components of v become zero (A = — oo), no information about the 




initial configuration may "percolate" to i = cx), and the asymptotic configura- 
tion is determined only by the random numbers used. This maximum Lyapunov 
exponent is also related to the synchronization properties of automata [17]. 

A preliminary numerical computation of A for our model is reported in Fig. 4. 
It can be noticed that the boundary j{R) of A > is not monotonic with R, 
reaching a maximum value for i? ~ 11. By comparisons with Fig. 2, one can see 
that the chaotic phase is included in the irregular one. 



4 Conclusions 

We have investigated a probabilistic cellular automaton model of an uniform so- 
ciety, with forcing local majority. The phase diagram of the model is extremely 
rich, showing a quiescent phase and several, active phases, with different dynam- 
ical behaviors. We have analyzed the properties and boundaries of these phases 
using direct numerical simulations, mean-field approximations and extending the 
notion of finite-distance Lyapunov exponent to probabilistic cellular automata. 
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